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Abstract. We consider electron-electron interaction effects in quantum point 
contacts on the first quantization plateau, taking into account all scattering processes. 
We compute the low-temperature linear and nonlinear conductance, shot noise, 
and thermopowcr, by perturbation theory and a self-consistent nonperturbative 
method. On the conductance plateau, the low-temperature corrections are solely 
due to momentum-nonconserving processes that change the relative number of lcft- 
and right-moving electrons. This leads to a suppression of the conductance for 
increasing temperature or voltage. The size of the suppression is estimated for a 
realistic saddle-point potential, and is largest in the beginning of the conductance 
plateau. For large magnetic field, interaction effects are strongly suppressed by 
the Pauli principle, and hence the first spin-split conductance plateau has a much 
weaker interaction correction. For the nonperturbative calculations, we use a 
self-consistent nonequilibrium Green's function approach, which suggests that the 
conductance saturates at elevated temperatures. These results are consistent with 
many experimental observations related to the so-called 0.7 anomaly. 
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1. Introduction 

Conductance quantization in a quantum point contact (QPC), first observed in 1988 p], 
constitutes a classic textbook effect of mesoscopic physics. On top of the integer 
conductance plateaus G = nGo (where Go = 2e 2 /h) observed as a function of gate 
voltage V g , many experiments have found a temperature-dependent suppression of the 
conductance appearing in the first half of the conductance plateau. This shoulder-like 
feature is seen at elevated temperature T (or finite voltage V) near the first quantized 
plateau[21 El H]> accompanied by a shot noise reduction^. It appears approximately 
around 0.7 Go, and has therefore been named the "0.7 anomaly" [21 El IH El El [3, E]- 
Despite the conceptual simplicity of a QPC and the fact that the 0.7 anomaly has been 
observed in a variety of material systems by different groups over more than a decade, 
still no generally accepted microscopic theory exists, apart from an overall consensus 
that one is dealing with some spin-related effect. 

While phenomenological models [9], assuming the existence of a density- dependent 
spin gap, can provide rather good fits to experimental data, the presumed static 
spin polarization due to interactions within the local QPC region is not expected 
in the presence of unpolarized bulk reservoirs. Along this line of thinking, it was 
recently conjectured that spin symmetry-broken mean-field or density functional theory 
calculations are unable to recover the correct T dependence of the conductance \TU\ [HI 
IT2] . A number of microscopic theories assume the existence of a quasi-bound state in 
the QPC region, leading to a Kondo-type scenario, as encountered in transport through 
interacting quantum dots[13l [H]. Such a quasi-bound state was indeed found in spin 
density functional theory (SDFT) calculations [13], and models based on this picture 
appear to reproduce several essential observations related to the 0.7 anomaly. However, 
different SDFT works also reached different conclusions [111 IT5] . Further proposals 
involve phonon effects[16]. Several publications have suggested that electron-electron (e- 
e) interactions alone may already result in a reduced conductance in a QPC at elevated 
temperatures, without the need for additional assumptions of spin polarization or a 
localized state pH HH [19l [201 EH E21 [23]. 

1.1. Main ideas and results 

Motivated by this body of experimental and theoretical work, we reconsider the role of e- 
e interactions for electronic transport properties in QPCs, starting from the assumption 
that in the low-temperature limit, a QPC is well described by a single-particle saddle- 
point potential. We then include e-e scattering, and, in particular, all momentum- 
nonconserving processes, where the number of left- and right-moving electrons does 
not have to be conserved in the scattering process, see figure [1]). These processes are 
not allowed in infinitely long translationally invariant (single-mode) quantum wires. 
However, the lack of translational invariance for a QPC connected to reservoirs permits 
such processes here. In other words, momentum is not a good quantum number 
for a QPC, and therefore interactions processes can violate momentum conservation. 
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Figure 1. Illustration of the two-electron momentum- nonconserving scattering 
processes that give rise to a correction to the transport properties at the beginning of 
the first plateau. The full (black) lines represent incoming electrons, while the dashed 
(red) lines are the outgoing electrons. The thick (blue) lines define the edge of the 
QPC. Only scattering between different spins is present to leading order in T/Tp due 
to the Pauli principle. 



Momentum-nonconserving processes are most relevant in the low-density regime, where 
the Fermi wavelength is comparable to the QPC's length, which is set by the curvature 
of the saddle-point potential and/or the distance to the gate electrodes [21]. Indeed, 
our quantitative analysis of the matrix elements for these processes (see below) shows 
that the effect of momentum-nonconserving scattering can be substantial, and implies 
that the conductance is significantly reduced at elevated temperatures, where the phase 
space for inelastic scattering is increased. We find that the breaking of translational 
invariance, and hence the backscattering rate, is most dramatic near the onset of the 
plateau, and then gradually decreases for larger electron density in the QPC. 

We start from the assumption that the low-temperature limit of a QPC at the 
first quantized plateau is well described by a Fermi-liquid picture with propagating 
single-particle states[25]. Throughout the paper, we consider only a single transverse 
channel being transmitted. Without interactions and at low temperatures, T <^ Tp 
(where e l F = ksT F is the local Fermi energy, see (fT8l) below), the conductance is given 
by the standard Landauer-Buttiker[26| formula, = GqTq(Ef)V , where V is the 
voltage difference across the contact and Go = 2e 2 /h. In this paper, we are interested in 
the properties at the plateau, i.e. when the zero-temperature transmission probability 
is close to one, 7o — 1. At zero temperature, electrons do not experience inelastic 
scattering from either phonons or other electrons. However, as temperature increases, 
phase space also increases for such scattering events. The effect of phonon scattering 
has previously been studied by Seelig and Matveev[16]. Here, we address the question 
of inelastic e-e scattering, which can cause electron backscattering. For example, two 
incoming electrons from the high-bias side can interact in the contact and scatter, so 
that one electron is backscattered while the other is transmitted. These processes are 
later denoted as fe(l) and shown in figure [Q Also shown is a backscattering process 
for two electrons |18j. which we denote as b(2). Both processes will lead to a current 
reduction. 

The matrix elements for the momentum-nonconserving processes can, to leading 
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Figure 2. Illustration of the effective ID potential of the QPC and the transmitting 
wavefunction ipE,ri- The wavefunction shows an enhanced weight in the contact region, 
which loosely speaking reflects the fact that electrons spend more time there. This 
so-called semiclassical "slowing down" is the main reason for the finite values of the 
backscattcring processes shown in figurc[lj cf. also [21], and will be discussed in section 
1231 



order in the single-particle reflection 1Z = 1 — %, be calculated using the fully 
transmitting wavefunctions. To good approximation, these can be described as WKB 
states in an effective ID potential Vq(x), which is a combination of the potential barrier 
and the confinement barrier [25], see section [2] and figure [2] The WKB states at energy 
E have the form 

VEtt ,(x) = y^f^exp (ir, J q dx'p{x')/^ , (1) 

where the local momentum is p(x) = y^2m[E — Vo(ar)], m denotes the effective mass, 
and rj = ± is the propagation direction. 

For a contact with T ~ 1, we may estimate the leading-order e-e interaction 
correction to the current by using Fermi's golden rule. To that end, we compute the 
rates for the two scattering events shown in figure [TJ and add them up, with a weight 
factor keeping track of the respective contribution to the current. This gives the e-e 
interaction correction to the current in the form (e > 0) 

I {2) = f E l^ 12 ' l'2')\ 2rh ' + riV ~ rh ~ ri2 (2) 

12,1'2' 

x mn 2 (l - ni')(l - n 2 ')5(E 1 + E 2 - E v - E 2 >), 

where we use the short-hand notation 1 = {E\, rji, o~\}, with a being the spin index, and 
sums run over the quantum numbers of the scattering states, 

£= ( 3 ) 

1 J cti=±t;i=± 

In ([2]), the occupation factors n\ are given by Fermi-Dirac distribution functions defined 
by the respective reservoirs, 

Ul = ^ ^ = e {E^ HR) )/k B T L{R) + 1 ' ( 4 ) 
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where one should use fii, for right-movers (r]i = +1) and /_t# for left-movers (r/i = — 1), 
and we have allowed for different temperatures in the two leads. Furthermore, in (T5]), 
the factor (771/ + 772' —rji — 1)2) accounts for the change of the relative number of right- 
and left-moving particles from initial to final states. Equation ([2]) is also found from 
rigorous perturbation theory in the Keldysh formalism, see section [31 where the matrix 
elements ^(12, 1'2') are specified in (PHI) . As expected, they contain both a direct and 
an exchange term, which is important for the behavior in a large spin-splitting magnetic 
field. The magnitude of the relevant coupling strengths is carefully estimated in section 



Based on (T5]), for To = 1, we obtain the linear conductance at low temperatures in 
the form 



where a realistic estimate for a typical GaAs QPC indicates that Ab ~ 1 at the beginning 
of the plateau. The result ([5]) has also been reported in [21] , where the QPC was 
modeled using a kinetic equation. Their prefactor Ab is proportional to an unknown 
"relaxation time in the leads", a quantity that does not appear in our theory. Instead 
Ab is directly connected to the inelastic e-e interaction processes. We find that Ab 
rapidly decreases when moving along the plateau, see also figure El The dimensionless 
coefficient A b includes the effects of both the b(l) and the b{2) backscattering processes. 
It is also interesting to consider the interaction correction at spin-split plateaus, i.e. 
for large magnetic fields. Here the energy dependence of the matrix element V^(12, 1/2') 
becomes important for equal-spin scattering, since exactly at the Fermi energy the direct 
and exchange terms cancel for equal spins only (see section 13721) . Hence, e-e scattering 
corrections in the spin-polarized case are much smaller. In fact, the leading contribution 
turns out to be of order (T/T F ) A , and the prefactor is smaller than Ab. Physically, this 
can be understood in terms of the Pauli principle. This qualitatively agrees with the 
experimental observation that almost no conductance suppression at the e 2 /h plateau 
occurs. 

Below we also provide results for the interaction corrections to other experimentally 
relevant quantities, such as shot noise, the thermopower, and the nonlinear conductance. 

The effect of e-e backscattering is thus most important at the beginning of the 
2e 2 /h plateau. On the other hand, for elevated temperatures, it leads to a breakdown of 
perturbation theory. Therefore, a crossover to a different type of behavior must appear. 
From equation ([5]), the temperature scale for this crossover is expected to be 



Contrary to the usual situation encountered in mesoscopic physics, the nontrivial 
question to be answered thus concerns the high-temperature limit (but still T <C T F ). 
To investigate this question, we have studied a local interaction model by means of a 
nonequilibrium formalism, employing a self-consistent Born approximation (SCBA) for 
the self-energy. In this case, our numerical results suggest that the linear conductance 



E3J 



G/G = 1 + G^/G = 1 - A b (irT/T F ) 2 , 



(5) 




(6) 
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approaches a saturation value g s = G/G of order g s ps 1/2 at high temperatures. The 
saturation is physically due to the fact that e-e scattering processes within the QPC 
region fully equilibrate outgoing electrons, which also suggests that g s is non-universal 
and thus depends on the detailed form of the various e-e couplings. We mention in 
passing that a high-T saturation of the conductance has also been reported for long ID 
wires H3 [27], 



1.2. Structure of the paper 

The structure of this article is as follows. In section [21 we define the model of the 
QPC, and provide estimates for the parameters involved. In section [31 this model is 
treated by lowest-order perturbation theory, and the interaction corrections to linear 
and nonlinear conductance, shot noise, and thermopower are computed. A simplified 
version of the QPC model with a local e-e interaction potential is then considered within 
a nonequilibrium formalism in section leading to a self-consistent numerical approach. 
This allows us to go beyond lowest-order perturbation theory, albeit in an approximate 
fashion. We briefly conclude in section [51 Details of the calculations can be found in 
two appendices. In intermediate steps, we sometimes set h — 1. 



2. Model and estimates 



2.1. Model 

We consider a two-dimensional electron gas (2DEG) with a single-particle potential U(x.) 
forming the QPC, where x = (x, y). Close to the middle of the constriction, x = (0, 0), 
the potential is assumed to be described by a saddle-point potential [28[ [29] 

1 1 

U (x) = U - -mu 2 x x 2 + -mu; 2 y 2 , (7) 

where x (y) is along (perpendicular to) the transport direction of the QPC. The QPC 
is thus characterized by the frequencies u x and ou y , or, equivalently, by the length scales 
£ y = \Jh/mujy and £ x = \J h/muj x . Transforming the 2D Schrodinger equation into ID, 
the ^-direction simply gives transversal ID subbands (modes) labeled by n — 0, 1, 2 . . .. 
Therefore the effective ID Schrodinger equation for mode n has a potential given by 

V n (x) = hu y (n + 1/2) + U - ^mu 2 x x 2 , (8) 
and the well-known transmission probability through the n mode is[28] 

T(E) = - (9) 

nV ; l + ex-p[-2ir{E - hu y (n + l) -U )/huj x Y w 

This transmission coefficient leads to the low-temperature conductance quantization in 
a QPC. 

In this paper, we focus on the situation when the Fermi energy is tuned such that 
only the lowest transversal state n = propagates through the constriction, and does so 
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with transmission close to unity, T (E E ) ~ 1. Therefore only the first (spin-degenerate) 
mode n = is included in the ID Hamiltonian H = H + Hj. The noninteracting part 
is 

H = £ J dx^{x) (-|^ + I'M, (10) 

and the Coulomb interaction gives 

Hi = ^Yl J dxdx '^t( x )^i'( x l W ( x ^ x ')^a'( x ')^cr( x ), (11) 

crcr' 

where ipl(x) (ip a (x)) is the creation (annihilation) operator for the n = mode. The 
effective (unscreened) ID interaction W(x, x') can be found by integration over the 
transverse eigenstates in the lowest mode, Xn=o(y) = exp(— y 2 /2l 2 ) / (n l l A ^JTy). This 
gives 

W(x, x') = ( - T =- M ( {X ~f )2 ) , (12) 

where k denotes the relative dielectric constant, M(x) = e x Ko(x), and Kq is the (zeroth- 
order) modified Bessel function of the second kind. This model of the interaction 
W(x, x') depends only on the difference x — x', and therefore appears not to break 
translational invariance. However, when including the scattering states, the effective 
interaction will break translational invariance, see the discussion in section [2731 below. 

For energies above the barrier top, the scattering states can be approximated by 
the WKB eigenstates ([TJ) , and we can thus expand the ID fermion operators as 

^M=Y1 dEl fE,ri(x)c Ejri:(7 , (13) 

T]=±l J 

where the energy integral is from the top of the single-particle potential to infinity, 
and c E a is a scattering-state annihilation operator with anticommutation relations 
{ c E,ri,ay c £v -q' a'} = ^v,v'^^'^( E ~~ ^')- The noninteracting Hamiltonian then effectively 
becomes 

H = Y,E 1 c\c 1 , (14) 
i 

where we again use the short-hand notation 1 = {Ei,r)i,o~x}, see fl3]). Similarly, the 
interacting part reads 



E WW,i2 cj, 4 c 2 c 1; (15) 

1,2,1', 2' 

with the matrix elements 

Wi'2',12 = 5 aiai ,5 a2(72 , J dx x dx 2 Ve 1/)J7i ,(xi)^ 2;j7?2 ,(x 2 ) 

x W(x 1 ,x 2 )cp EltVl (x 1 )<p E2tV2 (x 2 ). (16) 
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These matrix elements are discussed in section 12.31 For later purposes, it is also useful 
to introduce the corresponding matrix elements without the spin factors and taking the 
energies at the Fermi surface, 

W i% 2 ',r,im = J dx l0 lx 2 <p*E PtThl {xi)<P*E F , w M 

x W{x 1 ,x 2 )<p Epm {x l )<p EF!m {x 2 ). (17) 

Throughout the paper, we distinguish between the Fermi energy Ep in the equilibrium 
leads and the local Fermi energy ep in the QPC, 

4 = Ep — Uq — hu y = k B T F . (18) 

The single-particle thermal smearing of the conductance plateau occurs on the energy 
scale ep, and therefore we are here only interested in effects at lower energies, i.e. 

F- 



T <C Tf 



2.2. Hartree-Fock approximation 

In order to determine the saddle-point potential (j7j), one should in principle solve for the 
potential self-consistently, including the screening by surrounding gate electrodes and 
mean-field interaction effects. For the electrons in the constriction, the latter amount 
to the Hartree-Fock approximation, 

Hi jH f = ^ J ' dx J dx'VHF^ix,^)^^')^^), (19) 
where the self-consistent Hartree-Fock potential is 

V HF ,a(x, x') = S(x - x') I dx" W ( x > x")n a ,(x") 

- W(x,x')^t(x)M^)), (20) 

with n a {x) = (il)\(x)ip a (x)) . The self-consistent mean- field theory was discussed in [10] . 
showing that the mean-field result for the T-dependence of the conductance is markedly 
different from both the experimental observations and the finite-T corrections due to 
inelastic e-e processes studied in this paper. 

To illustrate this point, let us consider the Hartree approximation. Then the change 
of the local density n{x) with increasing temperature follows by using a Sommerfeld 
expansion. Noting that the chemical potential is set by the electrodes and can be 
assumed constant in this temperature range, we find to lowest order in temperature 



n{x) = n {x) 



1 



7T 2 / knT 



(21) 



24 \ep+ \muj^,x 2 . 

where no(x) is the T = density. First, we note that the density decreases with 
increasing temperature, because the ID density of states decreases with increasing 
energy. When (|2~T|) is inserted into the Hartree potential, i.e. into the first term in 
( j20|) . the temperature correction yields an additional contribution to the saddle-point 
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potential Vq(x). Expanding this result around the relevant barrier region x 
potential correction can be written as 

,2 



V H (x) = 5U 



A^ + 0((x/2E y r) 



0, the 



(22) 



Thereby the Hartree correction can be captured by a temperature-dependent 
renormalization of the barrier height Uq — > Uq + 5Uq{T) and of the curvature uj x — > 
\/l — A(T)cu x of the saddle potential in the transport direction. Explicitly, in terms of 
Meijer's G- function [30], we find 

3/2 



SUn 



with the energy scale 



W 



-W T 



fko n 



G 



3.2 
2,3 



U V e F 



1 1 



v / l87f47re K4 
Moreover, the dimensionless parameter A follows as 

2 



A 



W (k R T 



15 cu x hw x V Huj x 



(23) 



(24) 



(25) 



x 



G 



3.2 

2,3 



uj„e 



y c F 



huJ 2 x 



1 1 



-AG 



3.2 
2.3 



1 



The barrier is therefore lowered, 5Uq < 0, while the frequency u x is renormalized to 
smaller values, since (l25l) implies A(T) > 0. The quoted expressions are useful for 
7o — 1, where > huj x /2. Inspection of (Q then shows that the net effect of the Hartree 
correction is an enhancement of the transmission towards T = 1. As a consequence, a 
conductance plateau present already at zero temperature is then essentially not affected 
by the Hartree contribution as long as T <C Tp. A quantitative discussion of the Hartree 
corrections to the saddle-point potential can be found in section 12.31 below. 

Concerning the Fock part, a simple qualitative approximation is obtained if the 
pair potential is assumed to be a contact interaction, W(x,x') = Wq6(x — x'), see also 
[TO] . (This expression requires that the long-ranged tail of the e-e interaction potential 
is screened by nearby electrodes.) In that case, also the Fock contribution is local, and 
acts in precisely the same (but opposite) way as the Hartree potential, thereby partially 
cancelling the effects of the latter. In particular, it leads to the same temperature 
dependence of the corrections. We therefore conclude that the Hartree- Fock contribution 
leads to a very weak temperature dependence of the conductance, which even goes in 
opposite direction as compared to the effects of the inelastic e-e contributions. This 
temperature dependence of the mean-field result for the QPC conductance has been 
described in detail before [T0| [TTj. 

In the following, the weak T-dependence of the conductance under a single-particle 
picture will be neglected, and we focus on the effects of inelastic e-e scattering. Hence 
the interaction ( fTTl) is replaced by 

5Hj = Hj- Hj up. (26) 
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When performing diagrammatic expansions, the first-order contribution in 5Hj vanishes 
per definition, and the series starts with the second order in 5Hj, see sections |3] and HI 
From now on, the Hartree-Fock part is assumed to be contained in the single-particle 
potential Vq(x). 

2.3. WKB estimates for e-e matrix elements 

In order for our subsequent calculations to be meaningful, it is essential to first show 
that the relevant e-e backscattering strengths can be sufficiently large in experimentally 
relevant QPC geometries. In the following, we therefore provide such estimates and 
demonstrate that in the beginning of the first conductance quantization plateau, these 
couplings are indeed significant. Let us then consider an almost perfectly transmitting 
QPC and calculate the matrix elements that give rise to momentum-nonconserving 
scattering. For a simple estimate of the relevant e-e backscattering strengths, we now 
employ the WKB wavefunctions, see section 11.11 For larger reflection probability the 
WKB approximation is not adequate and a more sophisticated approach must be used. 

Because of the lower electron density in the QPC region, also denoted as 
semiclassical slowing downfIT\ . interactions are strongest near the constriction, see figure 
[21 The single-particle eigenstates in the potential Vq(x) are given by the WKB scattering 
states ([I]), which are a good approximation for %(Ef) ~ 1. The following estimates can 
thus be regarded as an expansion in 1Z = 1 — %, to leading order in 1Z « 1. In the 
WKB approximation, the matrix elements (fTTl) for e-e scattering at the Fermi surface 
are given by 



T/r/(o) - ( m V I A a < W ( X ' X ') 

Wwrvmm ~ { 2 nh) J p{x)p{x'\ 



p(x)p(x') 

x e i{m-vi>) So d vp(y)/ h e i{m-ri2') So d yp(v)/ h . (27) 

The semiclassical slowing down is reflected in the local density of states oc l/p(x), 
which is largest inside the QPC (x — 0). Here, p(x) = ^/2m[E F — V (x)} is the classical 
momentum, which we express using the ID saddle-point potential Vq(x) in (jEJ) as 

P(x) = -- — VT+WoY. (28) 
y y y 

The Fermi energy Ep enters here via the length scale a, 




fo-&V2) 



fUjJy 

This scale can also be related to the transmission Tq(Ef) by using ([9]), and hence to the 
position along the first conductance plateau. Inverting (Q, we find 

fY—^fi-V). ( 3o ) 

valid for 0.5 < % < 1. Next, consider the integral in the phase factor of ( 1271) . 

i(0 = \ J* dyv{y) = (3i) 
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tp 


rpt 
2 F 


O/iy 


%(E F ) 


fb(i) 


w b{1) 


A b (i) 


\hhJy 


2.3 K 


2.0 


0.985 


0.64 


720 eV" 1 


1.1 


fay 


5.2 K 


3.0 


1.0 


3.0 x 10~ 2 


34 eV' 1 


1.0 x 10" 2 




8.1 K 


3.7 


1.0 


1.0 x 10~ 3 


1.1 eV- 1 


1.5 x 10~ 5 



Table 1. Estimates for the interaction matrix elements Wb(i) and -A&m, see f|35[) . at 
three different positions on the first conductance plateau, i.e. versus the local Fermi 
energy e F . Parameters are: huj x = 0.3 meV, huj y = 0.9 meV, k = 10 and m = 0.067m e . 
See also figure [3J 



J- in 

2vr 



1-T 



where £ = x/a. This expresses (|27|) in a form suitable for numerical integration, 



{2nhxi x ) 2 Aire k£, 



~ fr)i' ,W , 



(32) 



where 



x M 



it - e? 



2 



(33) 



x cos[(?7i - rji'hiO + (m - V2')l(0] 

with M(x) defined after (1121) . The matrix elements for the 6(1,2) scattering processes 
in figure CD correspond to 



W, 



6(1) 



W 



(o) 



++>■ 



W, 



6(2) 



(o) 



while e-e forward-scattering processes are described by W&(o) = 
Wfi(i) enter the current through the dimensionless parameters 



W. 



(o) 



(34) 

The couplings 



,4 



16tt s 



6(i) 



0; 



12tt 2 



\ 2 



(35) 



We now use typical experimental parameters for a GaAs QPC:[3T] t\w x = 0.3 meV, 
hwy = 0.9 meV, k — 10 and m = 0.067m e . We have numerically integrated fb(i) at 
three different values of the local Fermi level ep, see table [U Clearly, W&m becomes two 
orders of magnitude smaller from the beginning of the conductance plateau to the end, 
while Aft(i) differs even by five orders of magnitude, see figure [31 

Let us then briefly discuss the magnitude of the Hartree corrections mentioned 
in section 12.21 Given the QPC parameters in figure El we obtain the energy scale 
W ~ 0.31 meV from (T241) . It is then possible to compare the Hartree conductance 
correction to the respective inelastic correction (jSJ). For concreteness, we take T = 0.7 K, 
corresponding to T* at the beginning of the plateau in figure [3J From 02 3|) . at the 
positions of the arrows in figure [31 we find 5Uq ~ —0.24 meV, —0.08 meV, and 
—0.05 meV (from left to right along the conductance plateau). Similarly, (I2"5j) gives 
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0.2 0.4 0.6 0.8 



Figure 3. The interaction matrix elements for backscattering of a single electron 
at various positions on the first quantized conductance plateau. The solid curve 
shows the T — conductance vs e F for a saddle-point potential with hw x = 0.3mcV 
and Hujy = 0.9meV. The decrease of the conductance in (j45|) is then orders of 
magnitude stronger in the beginning of the plateau. The backscattering of one 
electron (A^m) is significantly stronger than that of two (A^))- Neglecting A^ 2 ), the 
conductance corrections for the three values of ^Wi) arc: G^/Gq = — 2.1K~ 2 x T 2 , 
G( 2) /G = -3.6 x Kr 3 KT 2 x T 2 and G (2) /G Q = -2.2 x Kr 6 K -2 x T 2 , respectively. 
For T = T* , we have G^ = Go, and perturbation theory fails. This corresponds to 
T* = 0.7K, T* = 17K, and T* = 670K for the three stated values of A b(1) . Therefore, 
it is desirable to go beyond a perturbativc treatment. The quoted values are 

found numerically by integration of (|33|) . 

A pa 0.77,0.13, and 0.05, respectively. The Hartree corrections to the conductance are 
then captured by (jHJ), and imply an increase of the conductance. Since for all three cases 
in figure [31 the transmission at zero temperature is T > 0.985, the Hartree correction 
can at most enhance the conductance by a factor pa 1.015, which is a much smaller effect 
than the inelastic corrections discussed in this paper. Note that including the Fock term 
reduces the found Hartree contribution even further. 

Finally, screening by higher (closed) subbands and by nearby gate electrodes 
also leads to an interaction with broken translational invariance, where the effective 
interaction range is set by a combination of the pinch-off distance for closed subbands 
and the screening length of the electrodes. This would have two consequences, namely 
(i) an overall decrease of the e-e interaction strength, and thus of the Af,(i), and (ii) an 
additional source for broken translational invariance, which enhances the A^y The net 
result of including screening by the adjacent gates and 2DEG regions is therefore not 
straightforward, and a detailed electrostatic calculation is needed to gain a reliable 
understanding. Since semiclassical slowing down already results in very significant 
backscattering amplitudes, we focus on that in the present paper. 

3. Perturbation theory 

In this section, we compute the leading interaction corrections to the transport 
properties of an almost perfectly transmitting QPC. As discussed in section El we 
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treat the interaction Hamiltonian 5Hj by perturbation theory, yielding the current as 
J = + /( 2 ) + • - •. We here employ the Keldysh formulation to compute the correction 
1^ due to inelastic e-e scattering to second order in the interaction, and start with the 
definition of the current 

J = z e ^(lV(x)|l)^<(ll',tt), (36) 
ii' 

where Q < (ll',tt') = i(c\, (t')c 1 (t)) is the nonequilibrium lesser Green's function (GF), 
and J(x) = ^{d x — cQ. For T ~ 1, we can use the WKB states in (JI|) to calculate 
the current matrix element, 

(lV(x)|l) = ^A Wl(X)+W( " ) 



4tt^ y/pi{x)p v {x) 

x exp Q J^dylrjtptiy) - Vi'Pi'iv)]^ • (37) 

Here, we neglect terms proportional to d x p(x) to be consistent with the WKB 
approximation. In order to also allow for different temperatures in the two leads, we 
employ the nonequilibrium formalism in a slightly different way than usual[32J [55] . 
namely we include the different chemical potentials and temperatures in the initial 
density matrix p (H). (The left/right reservoirs inject electrons distributed by Fermi 
functions (J4]).) The interaction Hamiltonian 5Hi is treated as the perturbation that 
connects the two reservoirs. To evaluate the current, we thus write 

I = Tr\p(H )J{x,t)], (38) 

where J(x,t) is in the Heisenberg picture, J(x,t) = e lHt l h J{x)e~ %m l h . Here H = 
Hq + 5H T , and 

p(H ) = Z~ l e~( Ho > L ~^ LNL ^ kBTL ~( H °' R ~' J ' RNR)/kBTl1 , (39) 

where N^m and i?o,L/R are the number operators and noninteracting Hamiltonians, 
respectively, for electrons coming from the left / right lead. The noninteracting (retarded, 
advanced, lesser, and greater) GFs then follow in the form [32] 

g" la (i,uj) = (40) 

Qq(1,uj) = 2mni8{uj — Ei), 

Go(1,uj) = -2m[l-n 1 ]6(u-E 1 ), 

where n\ is given in (pEJ). The GF's are diagonal in spin space and identical for both spin 
directions. 

We now use a diagrammatic expansion of the interaction operator in terms of 
contour-ordered GFs[32]. Since the first-order diagrams vanish per definition, the 
leading diagrams are given by the second-order self-energy diagrams shown in figure HI 
The second-order correction to the lesser GF can be found by using standard Langreth 
rules|32J cf. (1731 . After some algebra, employing pUj) . we obtain 

./ ^ ^-^-V ^(i.")-w(i'.<-')i^ < (ii'.") 

ii' J 11 
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£(2) 



Figure 4. The second-order diagram included in The first-order and all 

other second-order diagrams belong to the self-consistent Hartree-Fock sector and are 
contained in the saddle-point potential, see section l^^l To avoid double-counting, they 
should therefore not be taken into account here. 

-£^(ii', cu)g<(i', W ) + g<(i, c)s( 2 ) a (ii', W ) } . (4i) 

One checks that this expression conserves current by differentiating with respect to 
x, and therefore we can set x — 0. Now there are two contributions to the current, 
corresponding to the real and imaginary parts of [E\ — Ey — iO 4 *]" 1 , where symmetry 
arguments show that the real part must vanish. For the imaginary part, we have a 
5-function enforcing E\ = Ey, which then, see (1371) . implies 771 = r]y. This leaves us 
with 

I {2) = -;|l> K 1 -*i)E (2)< (ll>£i) +n 1 ^ 2 >(ll,E 1 )] , 

where we use £( 2 ) r — Y>^ a = 'E^ > — E^ 2 ^. The lesser self-energy has two parts, the 
direct and the exchange contribution, corresponding to the two terms in figure HI In the 
time domain, this gives 

E®<(ll,tf) = -^^0 < (2,ttO^O > (3,ttO^O < (4,ttO(Wl3,24-W 3 i,24)W 2 4,l3.(42) 
234 

The expression for E^ 2 ^ follows by interchange of lesser and greater GFs in (|4"2I) . After 
a number of manipulations, we reduce this to 

j(2) _ _ ^ ^ |Wi2,l'2' - W21,l'2'| 2 g(?7l Vl> - V2') 

1,2,1', 2' 

x [n v n 2 '{l - ni)(l - n 2 ) - (1 - nv)(l - n 2 ')nin 2 ] 

x + £ 2 - E v - Ey). (43) 

In fact, this result recovers the Fermi golden rule expression (J2J), as may be seen by 
combining the two terms in (|43I) . We can now read off the matrix element entering (j2J), 

V (12, 1'2') = WV'2' - Wai.i/2' , (44) 

which involves the direct and the exchange interaction. 
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3.1. Conductance at low temperatures and voltages 

The perturbative interaction correction ( 1431) to lowest order in temperature T -C Tp 
and voltage eV <C s e F is determined only by the interaction matrix elements ( ITTj) taken 
at the Fermi surface. Some algebra then gives 



with the dimensionless coefficients Ab(i i2 ), see (j35l) . corresponding to 6(1) and 6(2) 
processes for electrons at the Fermi surface. We stress again that (l4"5l) holds only to 
lowest order in TZ Q = 1 — %, i.e. close to unity transmission at zero temperature. As 
described in section ITTT1 the process 6(2) corresponds to the simultaneous backscattering 
of two electrons, which has been discussed on a perturbative level in [18]. The process 
6(1), where a single electron is backscattered, has not been studied before. It is here 
important to note that for both processes, low-energy scattering is only possible between 
opposite spins: When the energy arguments in the matrix element (jUJ) are taken at the 
Fermi energy, the direct and exchange terms cancel each other for equal spins. The 
parameter A^i) is expected to dominate over Au 2 ), and has been estimated in section 
12.31 for a typical GaAs QPC We found that ^4&(i) changes from « 1 to ~ 10~ 5 from the 
beginning of the first quantized conductance plateau to the end, see figure [3] and table [TJ 
From ()45l) . we observe that interaction processes that do not change the number 
of left- and right-movers give no contribution to the current correction I^ 2 \ The 
same conclusion is reached for finite-length quantum wires using a Boltzmann equation 
approach [Ml 155]. In particular, forward-scattering processes do not contribute to l( 2 ', 
even though total momentum is not necessarily conserved in these processes either. Here, 
the important point is only whether the number of left- and right-moving electrons is 
conserved or not. However, at higher orders, important at elevated temperatures, all 
e-e interaction processes can come into play, see also Appendix A of [35j . 

3.2. Spin-polarized case 

Next we consider the effect of interactions for a completely spin-polarized QPC in a large 
magnetic field, i.e. on the e 2 /h plateau. For the spin-degenerate case, we found above 
that to leading order in T/Tp, only opposite spins interact due to the Pauli principle. 
As a consequence, the conductance correction is further suppressed for a single spin 
species, as we shall show now. In the spin-polarized case, we must take into account the 
energy dependence of the interaction matrix elements Wi 2t i> 2 >, which leads to a stronger 
power-law suppression in temperature. The current correction to second order in the 
interaction for a single spin species follows from (14"31 . and in linear response, we find 




(45) 




dE 1 dE 2 dE 1 ,dE 2 , (]V b{1)]] \ 2 + |H( 2 )||| 2 ) 



x nin 2 (l - n v ){\ - n 2 >) S(E 1 + E 2 - E v - E 2 >), 



(46) 
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where the rij are taken for fiR — [it — Ep. The same-spin interaction matrix elements 
are 

H(l,2)|| = W / Ei+T,E 2 +T>E 1 /±T,E 2 '-T ~ We x ^e 2 , 

i.e. a direct term minus an exchange term, where E\ and E 2 are interchanged. Here, the 
index Ey+ f — t] refers to the 6(1) [b(2)] process. To lowest order in temperature, 
the leading term is found by expanding the V&(i2)|i around the Fermi energy, since 
H(i,2)|| = when all energies are equal to Ep. The low-temperature interaction 
correction to the spin-polarized linear conductance = e 2 /h is obtained as 



h 



, (47) 



where the dimensionless parameter 



B bW = ^(k B T F T\(d El - d E2 )W Kl) ,\p F \ 2 {Ai 



with Wfc(i)-|- = We 1 +'[,e 2 +'[,e 1 ,+\,e 2 ,-'[ is estimated in |Appendix A| There, we find 
E>b(i) ~ 0.1 at the beginning of the quantized plateau (for %(Ep) = 0.985, corresponding 
to the first arrow in figure [3]). For the 6(2) process, the resulting power- law suppression 
is even of higher order, since 

(d El - d E2 )W bm \ EF \ 2 = 0, (49) 

where W 6 (2)t = Wei+\,e 2 +\,e xI -^e 2 ,-v To get (J49]), we use W 12 ,V2' = W 2 i,2'v, which in 
turn follows from W(x, x') = W(x', x). In conclusion, an extra factor (T/T F ) 2 suppresses 
the low-temperature correction at the spin-polarized e 2 /h plateau (i.e. the large- B case) 
as compared to the spin-degenerate 2e 2 /h plateau (i.e. the B = case). 

3.3. Noise 

Next we calculate consequences for another observable, namely nonequilibrium quantum 
noise. The zero-frequency shot noise follows from the (symmetrized) two-point 
correlation function of the current operator. Perturbation theory yields for the 
backscattering noise power in zero magnetic field 

S B (V, T) = 2e [2I bs{2) (V, T) coth{eV/k B T) 

+ J M1) (V, T) coth(eV/2A; B T)] , (50) 

where hs(i,2) are the current corrections due to Wb(i,2) quoted in (}4"5T) (defined positive 
for V > 0). This is nothing but the famous Schottky shot noise relation, encoding the 
charge of the backscattered particles. Equation ( l50l) predicts an additional factor of two 
for the b(2) contribution, because two electrons are backscattered in that event [T8] . 

The direct perturbative calculation of the noise to second order in SHj then yields 
the full noise power of the transmitted current in the compact form 

S T = S B + AG k B T - 8k B Td v I bsj (51) 
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where 1^ = Iwi) + /&s(2)- This perturbative result has also been obtained in a different 
context before[36]. In the limit V — > 0, one recovers the expected thermal Johnson- 
Nyquist noise 4:k B T[Go+G^ (T)] from ( JoTT) . with the interaction correction specified 
in Note that the last term in (j5T|) ensures that the correct thermal noise formula is 
obtained. 

Recent noise measurements on the first quantized plateau were compared to the 
corresponding single-particle picture[7], and a reduced noise power was observed at the 
conductance anomaly. Let us now connect ( f5Ti) to this discussion. For that comparison, 
according to [TJ, one subtracts the thermal noise and defines the excess noise as 

Sj = S T - 4G{V,T)k B T. (52) 

For a noninteracting system, this quantity is given by 

Sf p = 2G TZ[eV coth(eV/2k B T) - 2k B T] (53) 

to lowest order in the reflection coefficient TZ, see reference [TJ. In order to compare 
our approach to (1331 . we interpret Ibs/G V in a single-particle picture as an effective 
reflection probability TZ. Of course, this reflection is now mainly caused by interaction 
processes, and TZ is generally larger than the true noninteracting value TZq. Therefore, 
the difference between the true excess noise and its single-particle value is, in this 
framework, given by 

c _ cSP „y 



where 



h(x) = -8x + (vr 2 + x 2 )tanh(x/2). (55) 



From this expression, one observes that for eV < 6.507 k B T, regardless of ^6(1,2)) 
the measured noise should always be smaller than predicted by a single-particle 
analysis. This situation precisely corresponds to the parameter range of relevance for 
the experimental work of reference [7J, where eV < 5k B T. Our results for the noise 
are therefore consistent with the conclusions reached in reference |7]. Unfortunately, 
however, we cannot compare to the shot noise line profiles reported in Ref. 7, because 
our results are only valid to lowest order in TZq, and such line profiles would require a 
calculation function of 1Z . 

3.4- Thermopower 

Another experimental observable probing the enhanced phase space for e-e scattering 
at T ^ is the thermopower iS(T)[5], [37], for which perturbation theory for the spin- 
degenerate case predicts 

S(T) = ^ (A b(1) + A b{2) ) . (56) 
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Since the noninteracting thermopower is exponentially small [oc exp(— Tp/T)\ at 
the conductance plateau, the interaction correction completely determines the low- 
temperature thermopower [34] . The enhanced thermopower (as compared to the 
noninteracting one) is in qualitative agreement with experiments at the anomalous 
plateau [3 E]. 

4. Self-consistent nonperturbative scheme 

As discussed in section II. 1[ the intermediate temperature regime 

T* < T < T e F (57) 

needs to be understood beyond perturbation theory, see the caption of figure [3] for typical 
values of T*. Let us first discuss, on a qualitative level, what happens in limiting cases, 
where the QPC transport problem can be solved without invoking approximations. 
For instance, when A b t Qjl \ = but A b ( 2 ) ^ 0, one can establish correspondence to 
the problem of a single impurity in a Luttinger liquid with interaction parameter 
K = 2, which in turn is solvable [38] and leads to the high-temperature saturation 
value g s = G(T ^> T*)/Gq = of the linear conductance. Physically, this is quite clear: 
With strong 6(2) interactions, all particles are backscattered for increasing phase space, 
and the conductance approaches zero. Similarly, when A b ( lj2 ) = but ^ 0, i.e. 

when only forward scattering is present, the conductance is not affected at all, g s = % 
for 7o — 1. (However, for T < 1, forward scattering can increase the conductance, 
since such processes give an additive contribution to the usual elastic transmission.) 
Unfortunately, the most interesting special case, ^4&(o,2) = but A b ^ ^ 0, seems not 
described by an exactly solvable model. With only strong b(l) interactions present, all 
particles scatter in the QPC, but half of them are still transmitted. This suggests a 
reduction of the conductance by a factor of two, g s = 1/2. 

In order to make quantitative progress, we now consider a simplified model. It 
is given by the noninteracting part (fTOl) with Vq(x) = 0. Therefore, we now have the 
situation Ep = ep = Ef = ksTp. Moreover, we take the interacting part ( {Til with a 
local pair potential, see also [21], 

W(x,x') = W5(x)5(x'). (58) 

The local nature of the pair potential implies that all e-e interaction processes have 
the same amplitude, i.e. the amplitudes for backscattering of one and two electrons 
and the forward-scattering amplitudes are all equal, A b = A b ^ = A b ^ = A b ( 2 ). 
This is an artefact of the oversimplified local interaction (1581) . Nevertheless, a study 
of this problem is useful as it allows for an explicit calculation of the crossover from 
the T 2 corrections to the high-temperature conductance saturation. The model cannot 
be directly translated to the parameters relevant for a saddle-point QPC, because the 
dependence of the different backscattering matrix elements on the Fermi energy is not 
captured. Furthermore, because of the locality, the interaction potential (|58|) acts, in 
accordance with the Pauli principle, exclusively among opposite spins. Therefore, the 
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most important effect of a large 5-field, which consists of quenching the backscattering 
amplitudes, cf. section 13.21 is missed. We then only consider B = for the remainder 
of this section. For convenience, we express the interaction strength in terms of a 
dimensionless parameter, 

A = mW/2h 2 ir 3/2 . (59) 

Writing A in terms of Ab, we find[39] A = y/ZA^pK. On the first part of the plateau, we 
found Ab(\) ~ 1, see section [2TB1 with T*/Tp < 0.3. This allows for a meaningful study 
of the temperature range flBTl) with A < 1. 

In order to treat the local interaction (158]) . we start from the Dyson equation for 
the full Keldysh single-particle GF Q a (x,x']u)) for spin a, which is a 2 x 2 matrix in 
Keldysh space. Note that the interaction (1581) does not flip the spin, and therefore the 
Dyson equation is 

Q a {x, x'\ lo) = Q ;a{x, x'; lo) + Q - a {x, 0; u;)S ct (cj)^ ct (0, x'\ lo), (60) 

where the self-energy acts only at x = x' = 0. Moreover, both spins enter symmetrically, 
Q\ = Qii and we can suppress the spin index. In Appendix B we show that for the 
interaction ( 1581) and a parabolic dispersion, = k 2 /2m, the current can be expressed 
in terms of the local spectral function (at x — x' — 0) only, 

' = sEjf -/2M]£g' CD 

where the Fermi functions are defined in (j4j), and the local spectral function is 

A(cu) = i[g r {u) - g a {u)\ = -2 \mg r {u) (62) 

with Q{uj) = £7(0,0; uj). The noninteracting spectral function is 

A (lu) = 2txoI(uj) = (2m/ujf l2 6(uj), (63) 

where d(u) is the density of states, and 9 the Heaviside step function. Note that the 
real part of the noninteracting retarded GF Qq{oj) is nonzero at uj < for this dispersion 
relation, 

dk 1 



W J 2nuj-ky2m + t0- 

= -iid(u) (9(-u) +i9(u)). (64) 

In effect, the nonequilibrium current through the interacting QPC, ( 16Ti) . is fully 
expressed in terms of the local retarded GF only. 

So far, the given relations are exact, but to make progress, one needs to approximate 
the self-energy. As discussed in section 12.21 self-consistent Hartree-Fock diagrams are 
implicitly included in the single-particle potential ?7(x), so we are left with the higher- 
order diagrams representing inelastic e-e processes. Here we include the second-order 
self-energy in figure [51 made self-consistent by using the full GF. Note that the self- 
consistency ensures also current conservation |40] . This diagram amounts to treating the 
interactions effectively within the self-consistent Born approximation (SCBA). In the 
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Figure 5. The self-consistent second-order diagram for the spin direction a included 
in our SCBA scheme for the treatment of interactions. Double lines define full (self- 
consistent) GFs, as opposed to thin lines describing noninteracting GFs, see figure 2] 
Note that the exchange term exactly cancels the direct contribution for same spins, 
and therefore the second diagram in figure U does not appear here. 



perturbative calculations of section [31 the noninteracting GF enters instead, cf. figure HI 
In particular, the retarded component of the self-energy is 



where Qp denotes the local lesser/greater GF of the interacting system. The SCBA 
approach is similar in spirit to using a quantum Boltzmann equation approach with a 
self-consistent two-body collision integral, see Sec. 6.7.2 in [H] for a detailed discussion. 

4-1. Linear transport 

We first discuss the linear conductance, where the spectral function in (16TI) can be 
calculated in equilibrium. Thus, the lesser /greater GFs can be written in terms of the 
spectral function, 



where f°(uj) = \/[ e ^^F)/k B T + ±y Note that when eva i ua ti ng Q^(t) as needed in ( 1651) . 



negative frequencies have to be kept, cf. (J6l|) . Consequently, the self-energy and the 
interacting spectral function A(u) are generally nonzero at uj < due to the presence 
of interactions. However, the current only depends on the spectral functions at uj > 0, 
see (EH). 

To solve the self-consistency problem, we have employed an iterative procedure. 
Starting with the initial guess of the noninteracting GF, Q(w) = g (uj), see (j6l|) and 
( |6"6"|) . we compute £ r (cc>) from (1651) . which in turn defines a new retarded GF and a new 
guess for A{u) from the Dyson equation (|60|) . In the linear response regime, only the 
retarded part of (!60l) is needed, 




Q<(u) = %A{ U )f(u), G > (uj) = %A(u)(f(u) - 1) 



(66) 




(67) 



g r (u) = g r (cu) + g r (cu)z r (cu)g r (uj). 



(68) 
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With the solution to (1551) . we then compute a new estimate for X r (co>), and iterate the 
procedure until convergence has been reached. In the numerical implementation, it is 
essential to employ fast Fourier transformation routines to switch between frequency 
and time space, cf. (l67j) . This permits the fast evaluation of the self-energy ( )65|) in 
time space. For coupling strengths A < 0.8, this numerical scheme converges to a 
unique solution for the spectral functions A(oo), and can be implemented in an efficient 
manner. For larger A, however, several solutions may appear, and the approximation 
appears to be ambiguous. We therefore do not show results in this regime. Given the 
converged spectral function, we can compute the linear conductance G from (161]) by 
replacing — — > eV[— d u f (u)]. Thereby, G is numerically obtained as a function 
of the dimensionless parameters A and T/Tp. 

Let us first briefly discuss the lowest-order result for the retarded self-energy ( |65l) 
at T = 0. The result is obtained by first inserting ( 1661) with A(cu) —> A (uj), see ( 1631) . 
into ([67D. We find 



with Qo(t) following by overall sign change and 7 — > V, with the incomplete gamma 
functions j(a, x) and r(a,x)[5U]. For the chosen dispersion relation, we thereby get an 
UV-singular behavior in the self-energy E r (t) from (1651) . which manifests itself in the 
t — > behavior 



This divergence is also present at finite temperature, but affects only the real part of 
the self-energy. The perturbative results in section [3] are thus insensitive to it. However, 
when performing a nonperturbative calculation, a regularization of this UV divergence 
becomes necessary, and we chose a bandwidth cutoff of the order uj c « 3ep around the 
Fermi level. It is numerically convenient to employ a smooth cutoff function (e.g. a 
1/ cosh 2 filter), but the precise choice for the cutoff function does not appear to affect 
our results below. 

The numerical results for the linear conductance are shown in figure [6j First of all, 
we accurately recover the perturbation theory results at low T. At higher temperatures, 
a trend towards a conductance saturation could be conjectured from the numerical data. 
On a qualitative level, this saturation value could be expected from the following simple 
consideration. In the beginning of section H], we discussed that including only one of the 
backscattering processes of zero, one, or two electrons could lead to a high-temperature 
saturation value of g s = 2, 1, or 0, respectively. Since the point-like interaction gives the 
same amplitude to all three processes, it is tempting to simply take an average of the 
processes and conjecture a saturation value of e 2 /h for this model. 

Different functional forms for G(T) have been used to fit the experimental data 
for the 0.7 anomaly, including an activated T-dependence[3] or a phenomenological 
Kondo-type function We have tried to fit our numerical data to both types of 




(69) 
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Figure 6. Temperature dependence of the linear conductance G according to the self- 
consistent numerical approach described in section FOl for A = 0.3,0.6, and 0.8 (from 
top to bottom curve). The dashed curve gives the perturbative result for A = 0.8, dots 
refer to the numerical data, and solid curves are guides to the eye only. 




functions, slightly generalized to allow for different high-T saturation values. The 
activated behavior is 

® = l-(l-^)e- 1 !/ r 1 (71) 

where g Si \ characterizes the high-T saturation value, and T£ corresponds to the crossover 
temperature scale ([6]). The Kondo-like equation used for comparison with experimental 
data [I] is a scaled and shifted version of the universal scaling curve known from the 
Kondo problem. (However, this type of modification is not justified within the Kondo 
model, where G — > at high T.) We can get reasonable fits to both of these functions, 
see also [%T] . 

As a final remark on the numerical solution of the self-consistent approach in 
the linear regime, we mention that the thermopower (data not shown) exhibits a 
crossover from the S oc T 3 law at low T, see fl56l) . to a linear- in-T behavior at elevated 
temperatures. 



4-2. Nonlinear transport 

The iterative solution of the self-consistency problem is also possible out of equilibrium, 
and thereby allows to compute the nonequilibrium current flBTj) . The iteration now 
has to supplement fl65l) for Y7' by the corresponding equations for the greater /lesser 
components of the self-energy, 

/oo 
dte^g^(-t)gf(t)Gt(t). (72) 
-oo 

In addition, equation (1661) holds only in equilibrium, and has to be generalized. From 
the Dyson equation (loll , we obtain the retarded equation (Itjgj) and (for each a) the 
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Figure 7. Self-consistent numerical results for the voltage-dependence of the 
nonlinear conductance G = I/V for A = 0.3,0.6,0.8 at T = 0.02Tp (dotted curves, 
from top to bottom). For A = 0.8, also the comparison to the perturbative result 
(dashed curve) is shown. For the noninteracting case, G ~ Go for eV < 2sf (shown 
as black dash-dotted curve). 




relation 

g$(w) = gf{u) + £ »£ r H£%) (73) 

+ g r (u)^(u)g a (u) + gf(u)z a (u)g a (u), 

where the advanced components (g a , E a ) simply follow from the respective retarded 
ones by complex conjugation, since they are defined locally at x — x' — 0. For the 
iterative procedure, the initial values are now given by the lesser /greater noninteracting 
GFs 

g<(u) = i7Td(u)e(u)[f R (u) + f° L (u)}, (74) 

g>(u) = mdi^eiuj)^) + f° L (uj) - 2]. 

The iterative scheme can then be set up in a very similar manner, and yields a self- 
consistent solution for the local spectral function A(u), which in turn allows to compute 
the current from (16T|) . As a check, in the linear response regime of small bias voltage, 
we have reproduced the results of section 14.11 

Numerical results for the nonlinear conductance G(V,T) = I/V at a very low 
temperature are shown in figure [71 A clear decrease of the conductance with increasing 
V is observed, with a tendency towards saturation of the nonlinear conductance for 
large voltages. 

5. Discussion and conclusion 

In this paper, we have considered interaction effects in short QPCs. We have shown 
that taking into account moment um-nonconserving processes, we can qualitatively 
account for many of the experimentally observed behaviors of the linear and nonlinear 
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conductance (including their magnetic field dependencies), thermopower, and shot noise 
at the 0.7 anomaly. The gate- voltage dependence can also be qualitatively explained 
within the present scheme. The lack of translational invariance caused by a saddle- 
point potential is significant for realistic parameters often used to describe experimental 
realizations of QPCs. It is, however, also noteworthy that for longer quantum wires, 
where a saddle-point potential is not applicable anymore, similar backscattering effects 
can be caused by the lack of translational invariance at the ends of the quasi- ID wires. 

The self-consistent GF formalism for a simple point-like interaction model indicates 
that the conductance may approach a constant but non-universal value of order pa e 2 /h 
at high temperatures. The precise saturation value is expected to depend on the detailed 
e-e backscattering parameters, and hence on the detailed geometry of the QPC. The 
saturation should appear in the temperature dependence of the linear conductance. Our 
numerical data are also consistent with a saturation effect in the voltage dependence 
of the nonlinear conductance, but this issue requires further study. Physically, the 
saturation is caused by a relaxation of the incoming electron distribution functions due 
to the interactions within the QPC region. It should therefore be possible to develop 
a Boltzmann-type kinetic equation approach to describe (nonequilibrium) inelastic 
scattering in a QPC, see also [2TJ for attempts in this direction. However, such a 
development is not an entirely straightforward procedure and is outside the scope of 
this article. 

It is also instructive to discuss the point-like interaction model fl58|) in terms 
of an Anderson impurity model, see [JT] for details. By spatial discretization this 
Hamiltonian maps onto a ID tight-binding chain with hopping matrix elements t and on- 
site interaction U acting at one site (x = 0) only. In order to have perfect transmission 
at T = for the Anderson model also, the one-particle on-site energy must be canceled 
out so that the local energy is eo = Sp — ReS r (e^). We thus arrive at an Anderson- 
type impurity model similar to the one used in [13] to describe interactions in a QPC 
in the Kondo regime. However, we consider a rather different parameter regime where 
U is of the same order as the hybridization T and can be parametrically larger than 
the bandwidth D ~ \t\. Employing ([6]) with e e F ~ D, the interesting temperature 
range ( 1571) translates to D 2 /U < ksT <C D. While the Kondo model requires the 
formation of a local moment, this is not the case for the present approach. It would 
be interesting to explore this new parameter regime of the Anderson model by other 
non-perturbative schemes developed for the Anderson model. We believe that our 
theory offers a complementary approach to the 0.7 anomaly, which is not based on 
the assumption of a (quasi-)bound state in the QPC. In real samples, however, extrinsic 
effects (e.g. disorder) may be responsible for the presence of a bound state, and both 
mechanisms may then be relevant in parallel. 

To conclude, we hope that our work stimulates additional theoretical and 
experimental work on interaction effects in quantum point contacts. Future theoretical 
work should analyze the situation further away from the first conductance plateau, e.g. 
when T (E F ) is not close to unity |12"] . or the situation at higher quantization plateaus. 
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In the latter case, interaction effects due to semiclassical slowing are in principle also 
present on the higher conductance plateaus. However, in this case the other completely 
open channels screen the effective interaction, resulting in a presumably much smaller 
effect as compared to the first plateau. While we have repeatedly mentioned the 
experimental observations related to the 0.7 anomaly as the main motivation for our 
work, it is also clear that a realistic and quantitative description of experiments needs to 
consider a more refined modelling. This is especially important for the proper description 
of magnetic field effects away from the perturbative high magnetic field limit described 
in section 13.21 Moreover, on the theoretical side, it would be highly desirable to go 
beyond our perturbative calculation in a more realistic model. For instance, the gate- 
voltage dependence of the conductance at elevated temperatures - giving rise to the 
shoulder-like feature associated with the 0.7 anomaly - is not accessible to our theory 
at the moment, neither in the perturbative regime nor for the point-like interaction 
model in section IU where the gate- voltage dependence is difficult to incorporate. In any 
case, we hope that the model calculations discussed here can give new insights into the 
physics behind the observed anomalies in quantum point contacts. 
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Appendix A. Calculation of Bu X ) in the WKB approximation 



Here we estimate ( 1481) for the prefactor B^i) of the T 4 interaction correction to the spin- 
polarized conductance plateau e 2 /h, see ( 1471) . The large £>-field scattering amplitude 
involves a derivative of the WKB states with respect to energy, 

m 

,i I du- 
ll 



1 m i 
+ -rj 



2p 2 (x) 

When inserting this into ( 1481) . we obtain 

(d El -d E2 )W b{1)1 1 Ep 



dy 



p(y) 



2ixh/ 47renv27w£ 



x M 



K x — X 



r\2 
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2p 2 (x') 2p\x) 
Using the same parametrization as in section [21 this becomes 



h 

im 



dyp(y) 
dy 



p(y) 



(A.1) 



{d El - d E2 )W b{1)1 \ 



■ 9b(l), 



Electron- electron interaction effects in quantum point contacts 



26 



where 
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The value of Bu\\ is therefore given by 
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6(1) 



120tt 2 U 



x 8 / \ 4 



oo,, 



9b(i) 



Numerical integration for the same parameters as in section 12.31 (fku x = 0.3 meV, 
fkjy = 0.9 meV, k = 10, m — 0.067m e , and a/£ y = 2) then yields B b ^ « 0.1. 

Appendix B. Current for local interaction model 

The purpose of this appendix is to derive (ISTj) for the point-like interaction (I58I) . The 
derivation is a rewriting of the expectation value of the current operator, see (13"gi) and 

®, 

ft \ ^ 

2rm ' 



J(z,t) 



where x is arbitrary, ijj^{xt) is the creation operator in the Heisenberg picture, and the 
spin index is suppressed. In terms of the lesser GF Q < {xx' ,tt') = i(^(x't')xjj(xt)) , the 
average value is rewritten as 



(/) = -V r — I™ \d x > - d x ] g<(xx',uj). 
W 2m ^ 2vr L J CT v ' 7 



(B.1) 



Since the current is independent of x, we may take x = to evaluate it. To find 
d x Q < (xx' ,co) and d x 'Q < (xx', oo), the Dyson equation is used. The point-like interaction 
(I58p allows to simplify the Dyson equation from an integral equation to an algebraic 
one. Using the Langreth rules[32], we obtain for the lesser GF 

g < (xx', uo) = g$(xx', uo) + g r (xo, w)s r (oo, oo)g < (o x ', u) 

+ £ r (x0, cj)S<(00, uj)g a (0x', oo) 
+ £ <(x0, cu)S a (00, uj)g a (0x', oo). 

The noninteracting lesser and greater GFs g^(k, oo) include the Fermi functions of the 
leads. In ^-representation, they are 

g<(k,uj) = +2m5{oo - e k )f R/L (e k ), (B.2) 

g>(k,io) = -2m5(oo - e k )[l - f R/L {e k % 

where /^(sfe) is used for k < 0, and for k > 0. The Dyson equation for g < can 

also be written with the noninteracting GF £/o to the right of the self-energy, e.g. for g r , 



Electron- electron interaction effects in quantum point contacts 



27 



this looks like Q r (xx', uj) = Qq(xx', uj) + Q r (xO, a;)£ r (00, u)Qq(0x', uj). Differentiating the 
Dyson equation with Q on the left (right) of E with respect to x (x'), we get derivatives 
of the noninteracting GFs only, 

d x g< (oo, cu) = d x g<{oo,u) 
d x >g < (oo,u) = d x ,g$(oo,u) 

+ £ r (00, ^)£ r (00, uj) [3^(00, uj)}, 

where we use 

\imd x g r (xx',uj) = [ ^ - Zfc ^ =0, 



27T Co" — Ek + i0 + 

and likewise for the advanced GF d^^ixx' , uj), and for differentiation with respect to x'. 
Inserting this into IB. II) . and using g r {xx' , uj) = [g a (x'x, uj)}*, Y7(xx' ,uj) = [T, a (x' x , uj)}* , 
and d x Go(xx',uj) = —d x /go(xx',uj), we find 

f, /*0O I 

(/) = - / tR-^COO, w)] Re [1 + £T(00, w)57(00, w)] . (B.3) 

(J 

Furthermore, from ( 1B.2I) . one obtains 

/°° dk 
— g<(k,uj) ik 
-oo 21T 
777 

= e{oj)-[f R {uj)-f L {uj)]. 

Thereby we can express the current in terms of local GFs only, g(uj) = £7(00, uj). 
Combined with the Dyson equation g r = g$ + g^Y7g r at x — x' — 0, we obtain 

0(w)Re[l + = 0(uj)^t\ 

in terms of the local spectral function A(uj) = i[g r (uj) — g a (uj)]. After multiplication with 
— e < 0, we finally arrive at the electric current in (16T1) . Note that the current formula 
(|6ip can be seen as a continuous-space version of the Meir-Wingreen formula[43j for 
transport through an Anderson dot, see [JT] for a detailed discussion. 
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